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Abstract 

Background: Rhythmic clock widely occurs in biological systems which controls several aspects of cell physiology. 
For the different cell types, it is supplied with various rhythmic frequencies. How to synthesize a specific clock 
signal is a preliminary but a necessary step to further development of a biological computer in the future. 

Results: This paper presents a genetic sequential logic circuit with a clock pulse generator based on a synthesized 
genetic oscillator, which generates a consecutive clock signal whose frequency is an inverse integer multiple to that of 
the genetic oscillator. An analogous electronic waveform-shaping circuit is constructed by a series of genetic buffers to 
shape logic high/low levels of an oscillation input in a basic sinusoidal cycle and generate a pulse-width-modulated 
(PWM) output with various duty cycles. By controlling the threshold level of the genetic buffer, a genetic clock pulse 
signal with its frequency consistent to the genetic oscillator is synthesized. A synchronous genetic counter circuit based 
on the topology of the digital sequential logic circuit is triggered by the clock pulse to synthesize the clock signal with 
an inverse multiple frequency to the genetic oscillator. The function acts like a frequency divider in electronic circuits 
which plays a key role in the sequential logic circuit with specific operational frequency. 

Conclusions: A cascaded genetic logic circuit generating clock pulse signals is proposed. Based on analogous 
implement of digital sequential logic circuits, genetic sequential logic circuits can be constructed by the proposed 
approach to generate various clock signals from an oscillation signal. 

Keywords: Genetic oscillator, Clock signal, Logic circuit, Optimization 



Background 

Synthetic biology is an emerging interdisciplinary re- 
search field, which concentrates on understanding the 
behaviors of biological system from system-level as well 
as creating an artificial genetic circuit based on the prin- 
ciples of systems biology, mathematics and engineering 
[1-4]. Analogous to an electronic circuit, the synthetic 
genetic circuit also includes some standard biological com- 
ponents to assemble the biochemical process of living or- 
ganisms and achieve specific functionality. Based on a 
bottom-up approach, more complicated bio-computing 
modules can be expected to perform more complex func- 
tions via integrating a variety of biological devices, like 
very-large-scale integration circuits in electronics. By using 
mathematical models to capture the quantitative and quali- 
tative characteristics of biological systems, the customized 
genetic circuits with specific functions can be designed 
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from the system perspective [5-8]. For drug development 
and disease treatment, synthetic biology brings a useful 
and rapid direction through inserting the designed genetic 
circuits into the host cells to improve or modify the disease 
state of organisms. In addition, there are still potential ap- 
plications in biofuels, biotechnology, bioremediation, and 
bioenergy remained to be developed. 

Inspired by electronic circuits, several synthetic genetic 
circuits have recently been created, such as toggle switch, 
genetic oscillator, pulse generator, genetic counter, logic 
evaluator, sensor, filter, and cell-cell communicator. The 
former twos are based on protein-protein interaction with- 
out any external input to control their behaviors. Toggle 
switch applies two repressor genes repressing each other to 
cause bi-stable phenomenon, like as a memory device [9]. 
By cascading odd number of repressor genes in the cycle 
chain, a genetic oscillator can be synthesized to generate a 
stable oscillation signal in the protein response and applied 
in the control of dosage of drugs, or regarded as a syn- 
chronous mechanism for cell-cell communication [10-13]. 
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A pulse generator generates an instantaneous stimulat- 
ing signal and then resets to the original state by using 
time difference between the input and the corresponding 
delayed signals [14,15]. If the input signal is a periodic 
clock signal, then a clock pulse signal can be synthesized. 
Biosensor and filter are designed to detect the concentra- 
tions of specific molecular signal and range [8]. 

Boolean logic gate is an essential unit of a computer in 
digital logic circuits. To bring the insight of digital logic 
circuit design in electronic systems into biological sys- 
tems, the more complicated bio-computing processes 
can be easily constructed by combining a variety of gen- 
etic logic gates. The genetic logic gates constructed are 
based on different genetic transcriptional reactions to 
express various logic behaviors [16-22]. To use genetic 
components such as promoter, ribosomal binding site 
(RBS), repressor/activator genes and reporter gene, gen- 
etic logic gates with different logical operations have 
been assembled, such as NOT, Buffer, AND, OR, XOR. 
Through synchronous cascades of these genetic logic 
gates based on the topology of digital logic circuits, more 
complicated genetic logic circuits can be synthesized, such 
as multiplexer, half adder, combinational logic circuits, 
memory, and sequential logic circuits [23-30]. A genetic 
sequential logic circuit works with a counters, which is 
composed of some basic devices such as SR latch and flip- 
flop, has been developed in [25,29]. 

In biological systems, there are different rhythm frequen- 
cies depending on cell types. A 12-hour rhythm has been 
recently found in the mouse liver. For this reason, there are 
many engineered approaches proposed to synthesize the 
specific oscillation signals. In [31], the frequency-doubling 
oscillation can be constructed by using Fourier theory. A 
genetic circuit with multiple functions is designed to 
synthesize the oscillation signal with half original frequency 
[32]. Another aspect is to use regulated protein to control 
the transcription and degradation rates of target gene in an 
existing network structure [33-35]. A robust synthetic gen- 
etic circuit is designed based on control theory by regu- 
lating degradation rates of mRNAs and proteins in 
stochastic perturbational environments [33,34]. For cell-cell 
communication, synchronized genetic circuit designs are 
proposed to synchronize a population of oscillation signals 
[36,37]. To construct a promoter-RBS library from micro- 
array data and find suitable promoter-RBS components, a 
robust genetic circuit has been theoretically realized in the 
genetic systems by a systematic approach [38,39]. 

This paper proposes an artificial genetic sequential 
logic circuit with a function of frequency divider based 
on the periodic oscillation signal from a repressilator 
and analogous topology of the digital logic circuits in 
electronics. The proposed genetic sequential logic circuit 
is triggered by a clock pulse signal to generate a clock 
signal whose frequency is an inverse integer multiple to 



the genetic oscillator. Similar to an electronic waveform- 
shaping circuit, a genetic waveform-shaping circuit con- 
structed by several genetic Buffers in series is designed, 
which regulates time duration of logic high/low levels of 
an oscillation signal in the basic sinusoidal cycle and re- 
shapes the oscillation signal into a pulse-width-modulated 
(PWM) signal with different duty cycles by regulating the 
different threshold levels of the Buffer. The PWM signal 
can be regarded as a pulse signal with the frequency is co- 
herent to that of the genetic oscillator. The clock pulse 
signal is served as the rising or falling triggered edges of a 
clock signal with base frequency. In the digital logic the- 
ory, Karnaugh map is applied to determine the input 
signals of the rising or falling edge-triggered genetic JK 
flip-flops in each level [40]. A synchronous genetic 
counter circuit is triggered by the clock pulse signal to 
realize the genetic clock with its frequency is an inverse 
integer multiple to the genetic oscillator. 

For our proposed genetic pulse generator design, the 
periodic property of genetic oscillator is considered and 
the clock pulse signal is generated by utilizing the exist- 
ing synthetic genetic oscillator constructed by three re- 
pressor genes which repress each other in the closed 
loop. Different from the genetic counter circuit design 
[25,29], we introduce a generalized form based on the 
topology of digital logic circuits for synthesizing a clock 
signal with an inverse multiple of clock frequency to the 
genetic oscillator. The major advantage of the proposed 
approach is that it is easy to construct complex genetic 
sequential logic circuits via bottom-up approach with less 
computational time. Simulation results in silico show per- 
formance of the synthesizing clock pulse signal, and the 
clock signal with double, triple, quadruple basal periods 
while operating at the same genetic oscillator. 

Methods 

Dynamic model of synthetic genetic logic circuits 

By applying mathematical models to describe the bio- 
chemical reactions of genetic systems, a synthetic gen- 
etic circuit with a specific function can be synthesized 
from the system's perspective. 

Consider the dynamic model of the synthetic genetic 
logic circuit with L genes described by a class of nonlin- 
ear Hill differential equation [7] 

mi = aj \{u)-Ximi + a iiQ , ^ 
Pi= Pi m i-YiPi,i= 1,...,L 

where and p t denote, respectively, concentrations of 
mRNA and protein for the gene z, X t and y t are, respect- 
ively, the degradation rates of mRNA and protein, a t is 
the transcription rate of mRNA, fi t is the synthesis rate 
of protein, a iJQ is the basal production rate, is the 
promoter activity function which describes the nonlinear 
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transcriptional behavior and reflects the strength of the 
interaction between regulated protein and RNA poly- 
merase (RNAp), and u is the concentration of transcrip- 
tion factor (TF) which is produced from other gene(s) or 
inducer(s) to control the transcription rate of target 
genes. 

For a gene with an operator site which can bind a re- 
pressor or activator TF, the promoter activity functions 
are described as 



/not( w ) 



i+(#r 



(2) 



and 



/Buffer ( U ) 



or 



i + (f)" 



(3) 



where / NO t and^Buffer are promoter activity functions for 
logic NOT and Buffer [26,30], respectively, n is the Hill 
coefficient which denotes the binding cooperativity be- 
tween TF and the corresponding operator, and K is the 
Hill constant which are proportional to the lengths or 
affinities of the TF binding sites inserted into the pro- 
moter region of the target genes. For logic NOT gate, 
the input is a repressor and the gene produces a protein 
only in absence of the repressor; otherwise, the presence 
of the repressor obstructs the bound of RNAp and pro- 
moter. For genetic Buffer, the input is an activator which 
advances the bound of RNAp and promoter to produce 
protein. The frameworks for the two logic gates are illus- 
trated in Figures 1(a) and (b), respectively. 

For genes with two operator sites, which can bind two 
repressor TFs or activator TFs, the promoter activity 



functions are described in accordance with their logic 
functions as 
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where / AND , / OR , f XOR , f NAND and / NOR are, respectively, 
promoter activity functions of logic AND, OR, XOR, 
NAND and NOR gates, u x and u 2 are concentrations of re- 
pressor or activator TFs, K\ and K 2 are Hill constants for u x 
and u 2 , respectively, and n x and n 2 are the corresponding 




Figure 1 Expressions of a class of genetic logic gates, (a) NOT gate; (b) Buffer; (c) AND gate; (d) OR gate; (e) XOR gate; (f) NAND gate; and 
(g) NOR gate. 
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Hill coefficients. For logic AND, OR and XOR gates, the 
transcriptional behaviors are regulated by two activator TFs 
with different binding sites. Two repressor TFs control the 
genetic expressions of logic NAND and NOR gates. Their 
construction frameworks are shown in Figure 1(c)- (g). 

In [38,39], the promoter and RBS are considered as a 
promoter-RBS part to regulate the genetic expression 
because the half-life of mRNA is shorter than the corre- 
sponding protein has. One can rewrite (1) as 

Pi=pfi( u )-YiPi+Po,n i = !,...,£ (9) 
where 

Pi — i ' Pod — i 

Here, p t and p 0)i are new synthesis and basal produc- 
tion rates of the protein. The dynamic model of 2L dif- 
ferential equation (1) is reduced to the dynamic system 
with L differential equation (9). For real-world imple- 
mentation, fetching the corresponding promoter-RBS 
parts from the promoter-RBS library, the synthetic gen- 
etic circuit can be realized in the genetic systems. 

Synthetic genetic sequential logic circuits 

In digital logic circuits, the output of sequential logic cir- 
cuits depends not only on the present inputs but also on 
the past inputs. For synchronous sequential circuits, a 
clock signal is utilized as a metronome to coordinate ac- 
tions of circuits, which oscillates between high-level and 
low-level states. The circuits with triggered clock signals 
become active either in the rising edge, the falling edge, 
or in both of the rising and falling edges. For the se- 
quential logic circuit triggered at the rising edge of the 
clock signal, it becomes active when its clock pulse goes 
from low to high (0 to 1), and ignores high- to-low (1 to 0) 
transition. 

In genetic logic circuits, oscillation signal produced 
from a repressilator is not ideal as a clock for use in the 
kind of circuits relying on the change of rising or falling 
edge of the clock signal for state transition. Our pro- 
posed approach is to introduce the idea of a waveform- 
shaping circuit in electronics to genetic logic circuits, 
and reshape the synthesized genetic oscillation signal 
into a crisp clock signal or a PWM signal with different 
duty cycles. By regulating the size of duty cycle, the 
clock pulse can be generated with a rising edge or a fall- 
ing edge whose frequency is coherent to the oscillation 
frequency. To use the clock pulse, the designed genetic 
counter based on the topology of an electronic sequen- 
tial logic circuit is triggered to generate a clock signal 
with its frequency is inversely integer multiple to the 
genetic oscillation. 



Synthetic genetic oscillator 

Oscillation phenomenon in biological systems has been 
discovered at various levels of biological organization. Its 
practical function is to control the dosage of drugs or as 
a synchronous mechanism for cell-cell communication. 
The oscillation capability depends not only on the net- 
work topology but also on the system parameters. Cur- 
rently, the simplest synthetic genetic oscillator can be 
synthesized from a single gene repressing itself with a 
delayed negative feedback loop. An extension of the sim- 
plest oscillator, called a repressilator, consists of three 
genes {lad, tetR, cl) which represses each other in the 
cycle chain. The product of the first repressor gene, lad 
from E. coli, inhibits the transcription of the second 
repressor gene, tetR from the tetracycline- resistance 
transposon TnlO, whose protein product in turn inhibits 
the expression of the third repressor gene, cl from the A 
phage. Finally, cl inhibits lad expression, completing a 
negative feedback cycle [10]. The dynamic model of the 
repressilator can be described by 

Pi=PifNOT,i(P;)-YiPi ( 10 ) 

where p t and pj are concentrations of proteins for (/,;) = 
(lad, cl), {tetR, lacT) or {cl, tetR). For other design, the 
oscillation behavior can be generated by a number of re- 
pressor and activator genes in which the number of re- 
pressor genes must be odd. 

To design the genetic oscillator with desired oscilla- 
tions, one can realize a gene regulatory network to track 
a reference sinusoidal signal given by 

Vd = A sin(w 0 £ + <p) + y dfi (11) 

where y d is the oscillation signal with the desired ampli- 
tude A, basal frequency co 0 , phase cp and y dtQ is the base 
level to ensure nonnegative protein concentration. For 
more details regarding synthetic genetic oscillator design 
by optimization algorithms one is referred to [35]. 

Waveform-shaping circuit 

In electronics, a waveform- shaping circuit is designed to 
shape the input signal to the desired form according to 
an input and output (I/O) characteristic curve. For the 
oscillation input and the clock output, the I/O character- 
istic curve of the desired waveform-shaping circuit is 
displayed in Figure 2. A step function (dashed line) with 
a threshold level y T is used in electronics. For the input 
signal with its value larger than the threshold level, it is 
treated as "logic high". Otherwise, it is referred to "logic 
low". However, in biological systems, an ideal step function 
doesn't exist. A sigmoid function (solid line in Figure 2) 
might be used instead. From the I/O characteristic curve 
of a sigmoid function, there are two operational regions: 
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Figure 2 Idea of a genetic waveform shaping circuit. 

V J 



saturation and transition. The input signal in the saturation 
region can be cut-off and hold on the high level or the low 
level for approximation. In the transition region, the gain 
in the operation point y T must be more than (normalized) 
1 because it ensures that the input which is larger or less 
than the threshold level will be amplified or shrunk. By 
cascading the next sigmoid function, the oscillation input 
signal will gradually reach the saturation region and remain 
in the high or low level. 

According to this idea, a waveform-shaping circuit can 
be used to regulate the period of the logic high/low levels 
of an oscillation signal in a sinusoidal cycle and generate a 
PWM signal with different duty cycle defined by 

D = 100% (12) 
To 

where D is the duty cycle, T 0 is the basal period of oscilla- 
tion signal (11) with 2tt/co 0 and T on being the period of 
"logic high" in a basal period. For the PWM signals with dif- 
ferent duty cycles, the threshold is obtained by considering 

T 

y T =A s\n{o) 0 t + <p) + y dfi , t = t h ±-^ (13) 
with 

t h = — sm 1 (l)-^,t h e[0 T 0 ] (14) 

0)q COq 



To select the threshold level approaching to y dJQ + A, a 
clock pulse served as a rising triggered edge is generated 
and shown in Figure 3(a). For the clock pulse regarding 
as a falling triggered edge shown in Figure 3(b), one can 
choose the threshold level which is close to y d)0 -A. 
Similarly, a PWM signal with 50% duty cycle, i.e. the 
clock signal with its frequency is consistent to the gen- 
etic oscillator, is synthesized and shown in Figure 3(c) 
while selecting the base level of the signal y di0 . In other 
words, the PWM signals with different duty cycles can 
be synthesized from an oscillation signal via a waveform- 
shaping circuit in different threshold levels. 

Realizing a genetic waveform-shaping circuit 

In engineered genetic logic circuits, a genetic Buffer 
[30] is proposed to serve as a buffer between two cas- 
cade genetic circuits to enhance logic signal transfer. 
Its used here to aid the genetic waveform-shaping cir- 
cuit design: 

Pk = PhfBuf£er,k( u k,Kk, n k )-y k p k +p Qk , k = 1, ...,M 

(15) 

Its steady-state solution is easily obtained as 

PKss = ^fButterjMk,I<k,n k ) + k = 1, ...,M 

Yk Yk 

(16) 
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y T m ydfl~ A 



(c) 



Figure 3 Ideal PWM signals, (a) a rising triggered edge; (b) a falling triggered edge; and (c) 50% duty cycle. 



where p k is the output concentration of the /<th Buffer, p k>ss 
denotes its steady-state concentration, u k , K k and n k are, re- 
spectively, the input concentration, Hill constant, and Hill 
coefficient of the /<th Buffer and p k , y k an d Po,k are, respect- 
ively, synthesis, decay and basal rates. The second term of 
the right-hand side of (16) is the minimal level and p k /y k is 
the difference between the minimal and maximal levels. 
Output concentration of the genetic Buffer is the half max- 
imal output concentration when the input concentration 
equals K k and thus K k refers to the threshold level y T . 

In each stage, the corresponding inputs and the 
threshold levels are given by 



Pk-V 



k = l 
1 < k<M 



and 



f y T , k 

K k = I Pk-i +Po,*-i 



1 



2y, 



1 < k<M 



(17) 



(18) 



k-l 



In the first stage, the input signal is the oscillation sig- 
nal in (11) and the threshold level is chosen according to 
the desired duty cycle in (13). For the next stage, the 



input signal is the output concentration of the previous 
Buffer and the threshold level is the half maximal output 
level in the previous one. The topology of our proposed 
genetic waveform-shaping circuit is displayed in Figure 4. 
The oscillation signal from protein production of any 
gene of the repressilator activates the first gene in the 
genetic waveform-shaping circuit, whose production ac- 
tivates the next gene. Stage by stage, the oscillation can 
be reshaped to the crisp clock signal or PWM signal. 
However, the problem of slow convergence to the max- 
imal level is occurred for the larger threshold level K k . To 
resolve this problem, one can again cascade a Buffer with 
the design parameters of (16) in the last stage of genetic 
waveform-shaping circuit to compensate the output level. 
The gain at the operating point K k is obtained by 



du k 



u k =K k 



Pk^k 

*Y k Kk 



(19) 



where A k is the gain of the /<th Buffer. The gain is pro- 
portional to the Hill coefficient n k and the synthesis rate 
p k and is inversely proportional to the Hill constant K k 
and the decay rate y k at the operating point u k = K k . To 
ensure that the necessary condition of the gain at the 
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Figure 4 Topology of the designed genetic waveform-shaping circuit. 




operating point, K k should be exceeding 1. At first, one 
chooses the appropriate Hill constant for the desired 
synthesized PWM signal and then selects a suitable Hill 
coefficient n k) synthesis rate p k and decay rate y k satisfy- 
ing (19). From the system parameters in the previous 
stage, one proceeds to choose the appropriate system pa- 
rameters in the next stage satisfying (18) and (19). From 
[38,39], to realize the proposed genetic logic circuit in 
reality, one can find applicable promoter-RBS components 
from the constructed promoter-RBS library, whose I/O 
characteristic curves are capable of satisfying (18) and (19). 



Design of genetic frequency divider circuit 

Frequency divider in electronics is a device that gener- 
ates an output signal whose frequency is an inverse mul- 
tiple to that of the input signal. A sequential logic 
circuit, counter, is used to achieve this function, which is 
constructed by a series of flip-flops and triggered by the 
clock pulse to generate the clock signals with multi-fold 
basal period. Figure 5 illustrates an ideal clock signal 



while triggering at the rising edge of the clock signal 
with the desired basal period. 

Genetic JK flip-flop 

Genetic JK flip-flops based on the topology of digital 
logic circuits in electronics divide into the rising edge- 
triggered one and the falling edge-triggered one shown 
as in Figure 6. For the rising edge-triggered genetic JK 
flip-flop, its model is described by 

Pw = PwfAND(Pl<iPCLKi K Win W )-y w p w , 
Pv = PvfhNv{PjiPCLKi K Vi n v)-YvP V i 
PR = Pr/ ANd(PwiPqi K R^r)-YrPri 

Ps = PsfAND{PviPv K s,ns)-Y s Psi 
Pq= P(/nok^Pq^ k Q^q)-YqPq^ 
Pq= P-Qf^{Ps^PQ^ K ^ n - Q )-Y- Q P- Q 



(20) 



where Pclk is the concentration of clock pulse from low 
to high, p w , p v , p R , p s , p Q , and p^ denote, respectively, 
the protein concentrations of the genes W, V, R, S, Q 
and Q. The rising edge-triggered genetic JK flip-flop 
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(a) 



JUUL- 
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Pj P Q 

> JK-FF 

Pk Pq 



(b) 



uuir 

clock pulse 



Pj Pq 
> JK-FF 

Pk Pq 



Figure 6 A class of the clocked genetic JK flip-flops, (a) a rising edge-triggered one; and (b) a falling edge-triggered one. 



becomes active only when the clock pulse goes from 
low to high. There are four genetic AND gates and two 
NOR gates and the topology is displayed in Figure 7. 
The proteins p K and Pclk activate the transcription of 
the gene W. The proteins pj and Pclk activate the tran- 
scription of the gene V. The productions of the genes 
W and Q activate the transcription of the gene R and 
the productions of the genes V and Q activate the tran- 
scription of the gene S. The proteins p R and p^ inhibit 
the transcription of the gene Q and the proteins p s and 
Pq inhibit the transcription of the gene Q. 

For the falling edge-triggered genetic JK flip-flop, the 
model is described by 



Pw = PwfNAND(Pi<iPcLi<i K w,n w )-y w p w , 
Pv = PvfNAND(PjiPcLi<i K v,n v )-y v p Vl 
Pr = Pr/n and (p w i Pq i i n R)~ YrPr •> 
Ps = Psf nand (Pv>Pq> k s, n s ) -YsPsi 

Pq= pQfNAND(PSiPQi K Qi n Q)-Y Q PQi 
Pq~ Pq? NAND (PriPQ)Kq, Wq) ~YqPq 



(21) 



where Pclk is the concentration of clock pulse from high 
to low. This circuit is composed of six genetic NAND 
gates with its topological structure shown in Figure 8. 



Genetic counter 

To synthesize the clock signal with an inversely multiple 
of frequency of genetic oscillator, a synchronous genetic 
counter circuit can be used. The counter circuit in elec- 
tronics works on the rising or falling edge of the clock and 
count the number of clock pulses. Based on the feature, 
one first generates a series of clock pulses using our pro- 
posed genetic waveform-shaping circuit, and uses the 
clock pulse signal to trigger the genetic counter. Accord- 
ing to the Karnaugh map in the digital logic theory, the in- 
put signals of each genetic JK flip-flop and the topology of 
genetic counter circuit can be determined. To synthesize 
the clock signals with 2 s -fold basal period in which 6 is 
a positive integer, a synchronous genetic counter circuit 
with S number of rising edge-triggered genetic JK flip-flops 
is constructed by 



3 1 

?9 i 





> | Q Ui RBS | s 



Figure 7 Topology of the rising edge-triggered genetic JK flip-flop. 
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Figure 8 Topology of the falling edge-triggered genetic JK flip-flop. 



Pw y = Pw J and (Pk^Pclk, , K w, , ) -Yw^Pwx > 

P = Pv/AND (Ph'PCLK! > I< V l , «Vi ) "XV^V! > 

M = p*And , Pq x ? j ) ->v^i , 

Ps 1 = PS J AND (PV^PQS K Si > n S x ) -Ys l Ps 1 > 
^Qi = Pq/nor (Pi?x ^Qx > K Qi > W Qi ) "XQi^Qx ' 
>Qi = Pq/nOR fex , ^Qx >^Qx > "Qx ) -Xq 1 Pq 1 , 

Pvt. = PwJand(Pk s ^Pclk^ k w s ^ n Ws )-y Ws p w§ , 
Pv s = PvJand(P/ s ^Pclk^ k v s , n Vs )-y Vs p Vs , 
Pr s = PrJ 'and 0v 5 > Pq, > j n Rs ) -Yr s Pr s ' 

^S ff = Ps/aND (iV 5 > ^Q, > , «S 5 ) > 

^Q* = Pq/nor 0** '^Q, ' /C Q* ' n Qs)~YQ s PQs ' 
^Q, = Pq/nor (Ps, » Pq, » > "q, ) - Xq^q, > 

i*Gi = PGx/aND (p* 2 » ^Q 2 » 7 ^x , »Gi ) -Yg.Pg, , 

Pg,_ 2 = Pg s Jand{Pk s ^Pq s ^ k g s ^ n Gs _ 2 )-Y G8 _ 2 PG s ^ 

(22) 

with the input of each genetic JK flip-flop given by 



(23) 



where P C lk x * s tne dock pulse signal from low to high, 
PQi 1 PQi 1 PQs are ' res P ect ively, the clock signals with 
double, quadruple and 2^-fold basal periods. Figure 9 
shows the topology of the synchronous genetic counter 
for the clock signals with 2^-fold basal period. 



Ph 


=p K1 


= 1, 


Ph 


= Pk 2 


=P<h> 


Ph 


= Pk 3 


=p Gl , 


Ph 


= Pk s 


= Pg s . 



To synthesize the clock signal with triple basal period, 
the synchronous genetic counter with two rising edge- 
triggered genetic JK flip-flops and a falling edge-triggered 
genetic JK flip-flop is constructed by 

Pw,= PwJ 'and (Pk^Pclk, , K w, , n Wl ) -y Wl p Wl , 
Pv x = PvJand (PhiPcLid > K Vi > n Vxj-Yv x Pv x , 
Pr, = PrJ 'and (pwi > Pq, > K Ri > n Ri ) "X^i > 
^5x = Psx/and (p Vx 5 Pqx » 7 ^x , n Sl ) -r^x , 
^Qi = PqJnor fex ^Qx » /<c Qx » h Qi)~Yq 1 Pq 1 , 
^Qx = PQx/nor (Psx ,P Ql , ^q, , «q x ) -rq^Q, , 

i^^ 2 = Pw/and(Pk 2 ^Pclk^ k w 2 ^ W2 )-y W2 p W2 , 
Pv 2 = PvJand {Pj^Pclk, , K v 2 , ^k 2 ) -r^\/ 2 , 
^7? 2 = PrJ AND (Pw 2 1 Pq 2 i K Ri , ^ 2 ) -XtA , 
i^5 2 = Ps/aND {PV 2 1 PQ 2 i K S 2 , ^5 2 ) ~r5 2 ^5 2 j 

Pq 2 = PqJnor {Pr 2 i Pq 2 > /<C Q 2 ? W Q 2 ) -rq 2 ^Q 2 , 
^q 2 = PqJnor {Ps 2 iPq 2 > 7 ^q 2 > h Q 2 )-Yq 2 Pq 2 
Pw 3 = Pw Jnand (Pk 3 i Pclk 2 i K Wz ,n w ^) -y w ^Pw 3 , 
P k 3 = P v/nand (P/3 > ^ci/c 2 » >^ 3 ) -rv 3 ^^3 ' 

^7?3 = P7?/NAND (/V 3 ' ^Q 3 ' ^3 J W 7?3 ) -Xi? 3 ^3 ' 

= Ps/nand (^k 3 ' ^Q 3 ' /C5 3 j ) -Ys.Ps, i 
Pq 3 = PqjNAND(Ps^PQ^ K Q^ n Q,)-YQ,PQ^ 
Pq,= Pq /nand fe 3 > Pq 3 s ^q 3 , ^q 3 ) -ro 3J PQ 3 

^Gx = PGx/oR (^Q 2 , ^Q 3 , ^Gx , n Gx ) "XGi^Gx , 



with the input of each genetic JK flip-flop given by 

Pk x =Pk 2 = 1, 

P h =Pq 2 ,Pj 2 =P Ql , 

P h =Pq 2 ,Pks =Pq 2 



(24) 



(25) 
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where p C u< 1 * s the clock pulse signal from low to high, 
PciKi * s tn e clock pulse signal from high to low and p Gi 
is the clock signal with triple basal period. The topology 
of the synchronous genetic counter for the clock signal 
with triple basal period is displayed in Figure 10 and the 
corresponding ideal signals are shown in Figure 11. 

The above approach is generic, by an analogous way, 
one is able to determine the corresponding inputs of 
each genetic JK flip-flop based on the engineering digital 
logic theory [40] and cascade these basic flip-flips to re- 
semble other types of genetic counters with the desired 
operational frequency. 



where x is the regulation coefficient. The period of oscilla- 
tion increases when the regulation coefficient x decreases. 
Otherwise, the period decreases. Figure 12(a) shows the 
relationship between the basal period and the different 
regulation coefficient x< Oscillation with the basal period 
T 0 = 32 sec, the amplitudes A = 0.63, and the base level 
ya,o = 1.1731 is displayed in Figure 12(b) when x= 0.5. 

Synthesis of PWM signals 

Suppose protein concentration p cI is the desired oscillation 
input. To design the clock pulse signal with D = 10%, the 
designed genetic waveform-shaping circuit is described by 



Results 

To demonstrate the proposed synthetic genetic sequential 
logic circuit is effective to realize the function of frequency 
divider, the following numerical examples are illustrated 
to confirm the performance of the proposed method. 

Synthetic genetic oscillator 

Consider the dynamic model of the synthesized genetic os- 
cillator constructed by three repressive genes given by [34] 



Piad=X 0.685 



1 



0-233 PlacI Lpfc^O) = 0.7, 



PtetR =x( 0685 T^a 0-233/1 J, /^(O) = 1.2, 

p d = X ( 0.685 0.233^,^(0) = 1.7 



l + PtetR 



(26) 



Pi =fBuffer(PcI, 1 - 76 ^)-p l , 
P2 =/Buffer(/?lAM)-^ 
P3 =/BufferO2,0.5,4)-^3, 

Pi =/ B ufferfe O-M)-^ 
P5 =/Buffer0 4 A 5 5 4 )-^ 

p 6 = 1.0872/ Buffcr (p 5 ,0.5,4)- J p 6 



(27) 



where p 6 is the clock pulse signal with D = 10%. The cir- 
cuit has six genetic Buffers. In the first Buffer, the 
threshold level 1.76 is selected and the corresponding 
parameters including rate constants, Hill constants, and 
Hill coefficients and satisfying (18) and (19) are chosen 
for the second to the fifth Buffers. The last Buffer is to 
compensate the output of maximal level via choosing an 
appropriate rate constants in (16). In each stage, the I/O 
characteristic curves are shown in Figure 13(a) and the 
output concentrations are shown in Figure 13(b). 



(a) 



i 



> JK-FF 

r * ^ 



(c) 



Ph Pq* 
> JK - FF 



(d) 



> JK-FF 



(b) 




(f) 



Figure 10 Topology of the synchronous genetic counter for the clock signal with triple basal period, (a) Clock pulse signal from low to 
high; (b) Clock pulse signal from high to low; (c) Output signal p Q] of the first genetic JK flip-flop; (d) Output signal p Qi of the second genetic JK 
flip-flop; (e) Output signal p Qs of the third genetic JK flip-flop; and (f) Output signal of logic OR of (d) and (e). 
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(a) 



(b) 



(c) 



(d) 



(e) 



(0 



3T 0 

Figure 11 Ideal signals for synthesizing the clock signal with triple basal period, (a) Clock pulse signal from low to high; (b) Clock pulse 
signal from high to low; (c) Output signal p Q] of the first genetic JK flip-flop; (d) Output signal p Qi of the second genetic JK flip-flop; (e) Output 
signal p Qs of the third genetic JK flip-flop; and (f) Output signal of logic OR of (d) and (e). 



Similarly, the clock signal with D = 90% is generated 
by the following genetic waveform- shaping circuit 

Pi =/BufferO C /.°- 6761 . 4 )-A 5 
P2 =/BufferOlAM)-i^ 

P3 =J Bufferfe 0.5, 4)-/> 3 , ^ 

P4 = /Buffer(ft>°- 5 > 4 )-/>4> 
P5 = /Buffer(/>4>0.5,4)-/7 5 , 

^ 6 = 1.1083/ Buffe > 5 ,0.5,4)-^ 6 

where p 6 is the clock pulse signal with D = 90%. The first 
Buffer is designed with the threshold level 0.6761 and 
the second to the fifth Buffers choose the corresponding 
threshold level 0.5 satisfying (18). The Buffer in the last 



stage functions to compensate the output with the max- 
imal level. In each stage, the I/O characteristic curves 
are shown in Figure 14(a) and the output concentrations 
are shown in Figure 14(b). 

According the proposed approach, the oscillation 
signal can be rectified and filtered as a direct current 
(DC) output. The genetic waveform-shaping circuit 
with D = 100% is described by 

Pi =fBuffer(Pcn°- 5 ^)-Pn 

P2 =/Buffe>l,°-5,4)-^2, (?q \ 
^3=/Buffer(^2,0.5,4)-^, 1 > 

p 4 = 1.0851/ Buffer (p 3 ,0.5,4)-/7 4 , 
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where p 4 is the DC signal with D = 100%. There are four 
genetic Buffers. In the first Buffer, the threshold level 0.5 
is chosen and the second to fifth Buffers choose the 
threshold level 0.5. The Buffer in the last stage functions 
to compensate the output of maximal level. In each stage, 
the I/O characteristic curves are shown in Figure 15(a) 
and the output concentrations are shown in Figure 15(b). 

Synthesis of clock signals 

To design the clock signals with the double and the quad- 
ruple basal periods, we suppose 5 = 2, all synthesis and 
decay rates are 1, all Hill constants are 0.5, all Hill coeffi- 
cients are 4 in (22), p CLI<1 is the clock pulse generated from 
(27) and the high level DC signal is synthesized by (29). 



The concentration responses of the designed clock signals 
with double and quadruple basal periods are shown in 
Figure 16(a). State change of the genetic JK flip-flop occurs 
when the clock pulse goes from 0 to 1. For the clock signal 
with triple basal period, suppose that all variables remain 
the same to the above example, p C u< 1 an d Pclk 2 are > re ~ 
spectively, the clock pulses generated in (27) and (28), and 
the DC signal is synthesized by (29). Concentration re- 
sponse of the designed clock signal with the triple basal 
period is successfully generated, shown as in Figure 16(b). 

Discussion 

The goal of synthetic biology is to design genetic circuits 
with specific functions by using the approaches of 
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mathematics and engineering. Several genetic logic gates 
have recently been developed and experimentally real- 
ized [16-22]. In [16-18], simplified schemes of various 
genetic logic gates have been constructed by using basic 
biological components such as activator/repressor genes, 
reporter gene, promoter and RBS. For a gene with mul- 
tiple binding sites, the transcription rate of the gene can 
be regulated by the same number of regulated proteins. 
The modified protease can control the degradation rate 
of protein [19]. For logic NOT gate, the gene with a re- 
pressor regulated protein can generate the I/O charac- 
teristic curves with an inverse sigmoid function. A 
complementary form of NOT gate is "Buffer" [30] which 
is regulated by an activator protein to reshape the I/O 



characteristic curves for enhancing sharpness of a bio- 
logic response with a sigmoidal function. For logic AND 
gate with two binding sites, the genetic expression can 
be regulated by two activator proteins. In [22], the acti- 
vator proteins HrpR and HrpS control the promoter of 
the corresponding reporter gene to realize logic AND 
gate. By cascading a NOT gate in the outputs of AND/ 
OR gates, the NAND/NOR gates can be synthesized 
[20-22]. In [22], a reporter gene with different promoters 
(lac, BAD and lux) and different RBS parts (rbs30, rbs33, 
rbsH, etc) can determine the values of Hill constant and 
coefficient and achieve the different NOT gates with dif- 
ferent threshold levels. Similarly, one also can choose 
different promoter and RBS components with repressor 




input concent ran on (mM) time (sec) 



Figure 15 A PWM signal with D= 100%. (a) I/O characteristic curve in each stage; and (b) Concentration responses of the designed DC signal 
in each stage. 
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Figure 16 Response of the designed clock signals, (a) double and quadruple basal periods; and (b) triple basal period. 



or activator proteins to construct the corresponding sys- 
tem parameters for other logic gates. Through synchron- 
ous cascades of these genetic logic gates based on the 
topology of electronic circuits, the toggle switch and os- 
cillator have been characterized in real world by a class 
of Hill differential equations, and realized in Escherichia 
coli [9,10]. The previous papers have exhibited the possi- 
bility of realization of a class of fundamental biological de- 
vices in real world which form a basis for implementing 
more complicated combinational or sequential biological 
circuits. To realize the proposed genetic waveform- 
shaping circuit and genetic counter using a systematic ap- 
proach, we first can identify the different promoter-RBS 
parts to achieve the basic logic gates, and then assemble 
these logic gates with the topologies in Figures 4, 7, 8, 9, 
and 10. 

Reporter protein or green fluorescent protein (GFP) 
has been used in real-world experiments to reflect gen- 
etic expression [41]. By using a flow cytometer, the in- 
tensity of fluorescence of GFP is measurable. The system 
parameters in the dynamic model of genetic circuits can 
be identified from these measurement data using system 
identification methods [22]. In [38,39], a class of robust 
genetic circuits has been constructed by selecting the ap- 
plicable promoter-RBS component from a promoter-RBS 
library. Therefore, to realize our proposed genetic logic cir- 
cuits including the genetic waveform-shaping circuit and 
genetic counter, one first establishes the measurement de- 
vice which exhibits fluorescence concentrations of a series 
of a repressor or activator gene with different promoter- 
RBS components and TFs via fluorescence measurement 
[41] and then rebuild a promoter-RBS library with infor- 
mation of system parameters in terms of our mathematical 
model describing the behaviors of genetic logic gates by 
using the system identification methods [30,35]. To select 
the adequate promoter-RBS component generating the I/O 



characteristic curves of Buffers [30] and satisfying the de- 
signed conditions (18) and (19) from the rebuilt promoter- 
RBS library via an optimization algorithm (such as [35]), 
the clock pulse signals can be synthesized based on the 
proposed cascaded genetic logic circuit. Similarly, one can 
also choose the suitable promoter-RBS components to 
realize the various genetic logic gates and assemble these 
logic gates based the proposed topology. The designed gen- 
etic counters can be triggered by the synthesized clock 
pulse signals to generate the clock signals with multi-basal 
periods. 

In sequential logic circuits, the triggered signals could 
be rising or falling edge of a clock signal. One can use 
clock pulse signals to replace the rising or falling edge to 
trigger the sequential logic circuits such as our proposed 
genetic counters. The proposed genetic counters with a 
function of frequency divider based on the topology of 
electronic circuits can be used to generate clock signals 
with multi-basal period. In other words, one can utilize 
the proposed approach to generate clock signals whose 
periods is an integer multiple of 24-hour from a cell with 
circadian rhythm. 

Conclusions 

This paper has proposed a synthetic genetic sequential 
logic circuit as a frequency divider. The synthesized clock 
frequency is inversely multiple to that of the genetic oscil- 
lator which generates the fundamental sine wave. Through 
controlling threshold level of the genetic Buffer, the pro- 
posed waveform-shaping circuit regulates time duration of 
logic high/low levels in a basic sinusoidal cycle for an os- 
cillation input and generates ideal pulse signals with the 
coherent frequency to the genetic oscillator. Regarding the 
generated pulse signal as rising/falling edge of the clock 
signal with base frequency, a genetic synchronous counter 
circuit based on the topology of digital sequential logic 
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circuit is triggered by a pulse wave form to synthesize a 
clock signal with the inverse multiple frequency to the 
genetic oscillator. Experimental results in silico show the 
synthesizing genetic clock with double, triple, and quadru- 
ple basal periods while operating based on a genetic oscil- 
lator. By extending the proposed design principle, a class 
of multi-basal period clock signals can be generated in a 
straightforward manner. 
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